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Abstract. Power- law random banded unitary matrices (PRBUM), whose matrix elements decay in a power- 
law fashion, were recently proposed to model the critical statistics of the Floquet eigenstates of period- 
ically driven quantum systems. In this work, we numerically study in detail the statistical properties of 
PRBUM ensembles in the delocalization-localization transition regime. In particular, implications of the 
delocalization-localization transition for the fractal dimension of the eigenvectors, for the distribution func- 
tion of the eigenvector components, and for the nearest neighbor spacing statistics of the eigenphases are 
examined. On the one hand, our results further indicate that a PRBUM ensemble can serve as a unitary 
analog of the power-law random Hermitian matrix model for Anderson transition. On the other hand, 
some statistical features unseen before are found from PRBUM. For example, the dependence of the frac- 
tal dimension of the eigenvectors of PRBUM upon one ensemble parameter displays features that are 
quite different from that for the power-law random Hermitian matrix model. Furthermore, in the time- 
reversal symmetric case the nearest neighbor spacing distribution of PRBUM eigenphases is found to obey 
a semi-Poisson distribution for a broad range, but display an anomalous level repulsion in the absence of 
time-reversal symmetry. 



1 Introduction 

As a statistical description with few parameters, random 
matrix theory (RMT) [1] has found important applica- 
tions in a wide variety of topics, such as complex heavy 
nuclei [2], disordered electronic systems [3], large com- 
plex atoms [3] , stock market data analysis [5] , atmospheric 
data analysis 6; , analysis of human electroencephalogram 
(EEG) 7 , and spectra of complex networks [5], to name a 
few. Many interesting extensions of RMT have also been 
proposed to model different aspects of complex systems 
[9]. More relevant to this study, we note that RMT has 
been exploited extensively to predict the statistics of quan- 
tum systems with chaotic classical limits |10lllj . In these 
kind of systems, it is now almost confirmed that the near- 
est neighbor spacing distribution (NNSD) of energy levels 
follow the standard RMT prediction, namely, the Wigner- 
Dyson statistics [T^]- On the other hand, quantum sys- 
tems with regular classical limit, the same NNSD follows 
the Poisson statistics [13] . 

The (Anderson) delocalization-localization transition 
|14ll5j in disordered systems is a topic of enormous in- 
terest in condensed-matter physics because it is a metal- 
insulator transition induced by disorder. Recent advances 
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in cold-atom physics have also made it possible to directly 
observe this transition under controlled disorder |16ll7j . 
For weak disorder strength the eigenstates of the disor- 
dered system are essentially structureless and delocalizcd 
rather homogeneously over the whole Hilbert space. As a 
result the eigenstates can overlap very well with each other 
and level repulsion as a characteristic of the Wigner- Dyson 
statistics is anticipated. For strong disorder strength, the 
eigenstates are highly localized and their overlap is negli- 
gible. The energy levels are hence uncorrelated and level 
clustering as a feature of the Poisson statistics emerges. In- 
terestingly, in the vicinity of the delocalization-localization 
transition, the eigenstates display multifractal features [15] 
(thus neither delocalized nor localized) and the energy 
level statistics shows an interesting hybrid of the Wigner- 
Dyson statistics and the Poisson statistics [18]. To model 
such intermediate statistics for disordered systems under- 
going the delocalization-localization transition, the con- 
ventional RMT must be extended. Indeed, the power-law 
random banded matrix (PRBM) ensemble was proposed 
with success [15119] . A PRBM ensemble has two more pa- 
rameters than a conventional random matrix, and the ma- 
trix elements of a PRBM decay as a power-law function 
of the distance from the diagonal. 

The delocalization-localization transition and the as- 
sociated critical statistics may also occur in a quantum 
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system in the absence of any disorder, e.g., in a two- 
dimensional electron gas subject to a perpendicular mag- 
netic field |2QI21j or in driven quantum systems [22I23I17I24I 
For these cases, the well-established connection between 
classical chaos and the conventional RMT predictions, or 
between classical integrability and the Poisson statistics, 
breaks down. Of particular interest here is the statistics of 
the Floquet states of periodically driven systems. Indeed, 
critical eigenstate statistics and critical spectral statistics 
have been found in the kicked-Harper model [24], in an 
on-resonance double- kicked rotor model [55], and most 
recently, in a double- kicked top model [55], regardless of 
whether the underlying classical limit is regular or chaotic. 
These findings motivated us to propose in Ref. 57[ an ex- 
tension of Dyson's standard random unitary matrix en- 
semble [55J , such that the derealization- localization tran- 
sition in quantum mapping systems can be modeled with 
few parameters. Somewhat analogous to the Hermitian 
PRBM, we proposed in Ref. 27, a "power-law random 
banded unitary matrices (PRBUM)" ensemble, with the 
elements of each unitary matrix decaying as a power- 
law function of their distance from the diagonal. Remark- 
ably, as shown in Ref. |27j in terms of a fractal dimension 
analysis of the Floquet eigenstates, PRBUM around the 
localization-delocalization transition regime agrees with 
the statistics of an actual driven quantum system in the 
presence or absence of time-reversal symmetry, whereas 
the PRBM does not agree. This indicates that PRBUM 
is not a trivial extension of PRBM from Hermitian ma- 
trices to unitary matrices. This being the case, it be- 
comes necessary to investigate from many different an- 
gles the statistical properties of PRBUM. More specif- 
ically, we examine in this work the implications of the 
delocalization-localization transition for the fractal dimen- 
sion of PRBUM eigenvectors, for the distribution function 
of PRBUM eigenvector components, and for the NNSD 
of PRBUM eigenphases. Interesting findings include the 
dependence of the fractal dimension of PRBUM eigen- 
vectors upon ensemble parameters, comparison between 
PRBUM eigenvector component statistics with a normal 
distribution, transition of the NNSD from the Wigner- 
Dyson statistics to the Poisson statistics, and anomalous 
eigenphase level repulsion of PRBUM without time-reversal 
symmetry. Besides these driven quantum critical systems, 
the PRBUM ensemble has already found one application 
in a very interesting and important field of research. This 
ensemble has been found to be relevant to understanding 
sound propagation through underwater |29| . Our detailed 
results are hoped to motivate more interests in the poten- 
tial applications of PRBUM. 

This paper is organized as follows. In Sec. [5] we briefly 
introduce PRBUM first proposed in Ref. [57]. In Sec. [JJ 
we study the fractal dimension of PRBUM eigenvectors in 
the delocalization-localization transition regime as a func- 
tion of two ensemble parameters. Due to the finite size of 
the system we can afford to study computationally, it is 
not feasible to identify the exact transition point. For this 
reason we examine the transition point from different per- 
spectives in Sec. HI Statistical properties of the eigenvector 



components of PRBUM are examined in Sec. [5] In Sec. 
[6] we seek the signature of the delocalization-localization 
2&/2fisition in the NNSD of PRBUM eigenphases. Finally, 
we conclude this paper in Sec. [5] 

2 From Dyson's circular ensembles to 
PRBUM ensembles 

Dyson introduced three types of random unitary matrix 
(circular) ensembles to study statistical properties of com- 
plex quantum systems [28]. Dyson's circular ensembles 
have also proved to be useful to study periodically driven 
quantum systems whose classical limits are chaotic. In par- 
ticular, the circular unitary ensemble (CUE) models com- 
plex quantum systems without time-reversal symmetry. 
The CUE is an ensemble of (N x N) random unitary ma- 
trices distributed with the natural Haar measure on the 
unitary group U(iV). The Haar measure is analogous to a 
uniform distribution which implies that the probability of 
having any unitary matrix in the CUE is equal. The eigen- 
phases of CUE are hence uniformly distributed on the unit 
circle. The circular orthogonal ensemble (COE) is an en- 
semble of symmetric unitary matrices that model systems 
with time-reversal symmetry. A COE matrix U can be ex- 
pressed in terms of another unitary matrix W from CUE, 
i.e., U = W T W . In addition, sampling a matrix U from 
COE is equivalent to sampling a matrix W from CUE [50] , 
The third type of Dyson's circular ensemble is not related 
to this work and hence will not be discussed here. 

One convenient numerical algorithm for the generation 
of CUE and COE is Mezzadri's algorithm [31]. The first 
step in this algorithm is to sample matrices from an en- 
semble of N x N random matrices {Zq}, with both the 
real and the imaginary parts of the matrix elements be- 
ing Gaussian distributed random numbers (such a com- 
plex matrix ensemble is known as the Ginibre ensemble 
in the literature). As already detailed in Ref. [57], the key 
element in our algorithm for PRBUM generation is to re- 
place {Zq} by matrices drawn from the PRBM ensemble, 
with all other following steps being still the same as in 
the original Mezzadri's algorithm. The PRBM ensemble is 
composed of random Hermitian matrices whose matrix el- 
ements {Hij} are sampled from independently distributed 
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Fig. 1. (Color online) The parameter ao is plotted as a 
function of the parameter a for (a) PRBUM-COE and (b) 
PRBUM-CUE. Here b = 0.1 (•), 0.3 (A), and 0.5 (■). 
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Since PRBM has two parameters a and b (the case of 
a = 1.0 represents the critical Anderson transition point 
in a time- independent problem), the PRBUM generated 
from PRBM can also be characterized by a and b. Par- 
allel to the CUE and COE matrices generated from Mez- 
zadri's algorithm, our modified algorithm generates CUE 
and COE versions of PRBUM ensembles, which are called 
PRBUM-CUE and PRBUM-COE below. 

Other details about the generation of PRBUM can be 
found in Ref . (27) . Remarkably, the variance of the matrix 
elements of the generated unitary matrices does display a 
power- law feature, i.e., 
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In Fig. [TJ we present the parameter a>o of the PRBUM 
ensemble as a function of the parameter a of the initial 
PRBM ensemble. This result shows that the numerically 
found values of «o depend on the parameters a and b of 
the PRBM used in the first step of our algorithm. More- 
over, even with the same set of a and b, the values of 
a for PRBUM-CUE differ from those for PRBUM-COE. 
For this reason, we still prefer to use the parameters a and 
b to characterize our PRBUM. Figure [2] depicts the nor- 
malized density distribution of the eigenphases for both 
PRBUM-COE and PRBUM-CUE. Eigenphase 9 of a uni- 
tary operator U is defined in the following way: 



U\ct>x) = expWxMx) 



(3) 



where {|0a)} is an eigenstate of U. The distribution is 
found to be uniform as in the CUE and COE cases. Here 
we show result only for b = 0.1 with a = 0.5. For other 
values of a, we also find similar uniform density. This re- 
sult is important because it indicates that the PRBUM 
ensemble also satisfies Haar measure and our algorithm 
is very robust. As such, there is no need to carry out an 
unfolding procedure when we study the NNSD. 



3 Fractal dimension analysis of PRBUM 
eigenvectors 

The fractal statistics of PRBUM eigenvectors can be an- 
alyzed by the inverse participation ratio (IPR), 
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where A is the index for eigenstates \<f)\) and {|m)} are 
generic basis states. In our early study [27] . it was found 




Fig. 2. (Color online) Normalized distribution [denoted 
p(6)] of the eigenphases 9 for (a) PRBUM-COE and for 
(b) PRBUM-CUE. Here the parameters b = 0.1, a = 0.5. 
The matrix dimension N — 50 and the ensemble size is 
10000. Numerical results here indicate a uniform eigen- 
phase density. 

that for a = 1.0, the IPR scales anomalously (i.e., D 2 is 
fraction) with the Hilbert space dimension N, i.e., 
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where D% is just the fractal dimension of a particular 
eigenstate \4>\). Following the methodology in studies of 
the PRBM ensemble [32 33 , it was further observed in 
Ref. [27| that for sufficiently large N, the distribution 

of ln(Pj A ^) for a = 1.0 displays evidence of being scale- 
invariant. This scale-invariance, which also suggests self- 
similarity, is valid upto a certain scale. We define D 2 in the 
non-critical regime where, at the thermodynamic limit, 
the system should be either in delocalized (D 2 = 1.0) or 
localized (D 2 = 0.0) regime. However, due to the bounded 
scale, some part of the delocalized or localized regime 
(close to the critical point) may show self-similarity like 
at the critical point with different fractional value of D 2 . 

Basically the technique of determining D 2 is follow- 
ing: due to the scale-invariance, the distribution function 
of ln(i- > 2 A ' ) ) only shifts as N varies. As such the spectral 

average of ln(P 2 ), denoted by (ln(P2)), is linearly related 
to IniV. The slope of the (ln(P2)) vs In TV curve gives a 
single fractal dimension D 2 for the system. The scaling of 

the variance of P 2 (A) with N is also studied in Ref. [27] 
for a = 1.0, with features different from those observed in 
PRBM. We have also observed that conclusions drawn in 
Ref. 27 based on PRBUM for a = 1.0 are also valid if we 
apply our analysis to a different random unitary matrix 
model that also aims to describe fractal and intermediate 
statistics in quantum mapping systems [34135] . 

The purpose of this section is twofold. First, we present 
further numerical evidence that like PRBM, the case of 
a = 1.0 is at least (if not exactly) close to the critical 
transition point of PRBUM. We do so by examining the 
transition of the fractal dimension D 2 as a function of a 
for various values of b. Second, we examine the depen- 
dence of the fractal dimension D 2 upon b. Note that since 
we are not in the thermodynamic limit, so the D 2 values 
only represent some self-similar properties within a certain 
scale. 
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In Fig. |31 we show the dependence of the numerically 
obtained D2 as a function of a for three different values 
of 6, i.e., b = 0.1,0.3 and 0.5. Similar to Ref. [37], we 
have considered matrix sizes ranging from a few hundreds 
to a few thousands to extract the value of Z?2- Left and 
right panels of Fig. I are for PRBUM-COE and PRBUM- 
CUE, respectively. For both cases, as the value of a scans 
through unity, we have the change from D2 ~ 1.0 to 
£>2 ~ 0.0, which reflects the transition from extended 
states to localized states. This is always the case despite 
the big change in the values b (see discussion below). In 
this sense, our PRBUM ensemble may be regarded as an- 
other unitary model of the Anderson transition |36j . How- 
ever, the delocalization-localization transition is not very 
sharp; eigenstates are fractal (1.0 > D2 > 0.0) for a some- 
what wide range of values of a around a = 1.0. We will 
discuss more on this issue in the next section. 

In Fig.UJa) an d (b) we show how the fractal dimension 
D2 of PRBUM eigenstates for a = 1.0 depends on b. For 
the PRBUM-COE case, in b < 1 regime, D 2 is found to 
change linearly with 6, namely, D2 — 36 (dashed line in 
Fig. Ufa)). This feature is different from the parallel result 
from the PRBM ensemble at a = 1.0, which gives D2 = 2b 
[33 . One may wonder if this difference is simply due to the 
difference between b and 60. So in Fig.HJc), we also show 
how D2 depends on bo . Here we like to mention that, when 
b varies from 0.01 to 0.2, correspondingly bo varies from 0.2 
to 0.8. In this regime, we find D2 — 1.136o (dashed line 
in Fig. [He)), this expression is giving a proportionality 
constant far from two. For the PRBUM-CUE case, in the 
same regime, D2 still varies linearly with b (dashed line 
in Fig. Bib)); but the proportionality constant is not an 
integer (D2 — 2.686). On the other hand, D2 varies with 
60 as D 2 — 2.226o (dashed line in Fig. Efd)). By contrast, 
in the regime 6 6 [0.2, 1.0], D 2 approximately follows 6 as 
D 2 = a — c/b (solid line in Fig. H (a) and (b)) where a 
and c are constants. We also find from our numerical fit 
that, for the PRBUM-COE case a ~ 1.11 and c ~ 0.10; 
whereas for the PRBUM-CUE case, a ~ 1.17 and c ~ 
0.12. This behavior is qualitatively similar to the behavior 
observed in case of the PRBM ensemble. For PRBM in 




Fig. 3. (Color online) The fractal dimension D2 of the 
eigenvectors of PRBUM plotted as a function of a, for 
6 = 0.1(»), 0.3(A), and 0.5(B). For both (a) PRBUM- 
COE and (b) PRBUM-CUE, one may identify that the 
delocalization-localization transition occurs around a — 
1.0. 
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Fig. 4. (Color online) The fractal dimension D 2 of the 
critical eigenstates of (a) PRBUM-COE and (b) PRBUM- 
CUE as a function of 6. For the PRBUM-COE case, in 
the 6 <C 1 regime D2 varies linearly with 6 as D2 = 36 
(dashed line). In the regime 6 € [0.2,1.0], D2 approxi- 
mately follows D2 = a — c/b (solid line) where the param- 
eters a ~ 1.11 and c ~ 0.10. Finally, in the 6 > 1.0 regime 
D2 rapidly saturates at a value close to unity (dashed- 
dot line). For the PRBUM-CUE case, D 2 behaves qualita- 
tively the same as in the PRBUM-COE case. However, in 
the 6 1 regime, D2 does not exactly follow D2 = 36. In 
the 6 € [0.2, 1.0] regime, D2 again approximately follows 6 
as D2 = a — c/b with a ~ 1.17 and c ~ 0.12. The same D2 
is also plotted as a function of 60 for (c) PRBUM-COE and 
(d) PRBUM-CUE. There, for the PRBUM-COE case and 
for 6 G [0.01,0.2] and hence 6 € [0.2,0.8], D 2 ^ 1.136 
(dashed line). In the regime of intermediate 6 values, we 
cannot fit D2 vs 60 curve by the function D2 = a — c/bo- 
For even larger 60 values, similar fast saturation of D2 is 
observed. The PRBUM-CUE case behaves qualitatively 
the same. 



the same regime, D 2 = 1 — 1//3tt& where (3 = 1 and 2 
for the orthogonal and unitary symmetries, respectively 
[3"3"] . Note also that we cannot get a nice fitting with the 
same form for the D2 vs. 60 curve. Furthermore, in the 
regime of 6 > 1 , D 2 is found to rapidly approach a constant 
value 1 - e, where e ~ 0.013 for the PRBUM-COE case 
and e ~ 0.011 for the PRBUM-CUE case (dashed-dot line 
in Fig. IDJa) and (b)). Similar fast saturation of D 2 to a 
constant value is also observed in D2 vs. 60 curve. Results 
in Fig. 0] also suggest that tuning the 6 value from 6 = 0.1 
to 6 = 0.5 as in Fig. [3] already covers a significant range. 



4 Is a = 1.0 the transition point? 

In the previous section, we have shown that the PRBUM 
delocalization-localization transition takes place somewhere 
around a = 1.0. The exact transition point is not deter- 
mined. Numerically, similar results were also observed in 
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Fig. 5. (Color online) Ensemble average of the frac- 
tal dimension di is calculated using box counting pro- 
cedure as a function of the system size N for a — 
0.80,0.90,0.95,1.00,1-05,1.10 and 1.20 from top to bot- 
tom. Here we construct the ensemble by taking 1200 to 300 
matrices for different Hilbert space dimensions. The re- 
sults for PRBUM-COE and PRBUM-CUE are presented 
respectively in panel (a) and (b), with b = 0.1. 
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Fig. 7. (Color online) The slope ?y of g?2 vs. log(iV) curves 
(as presented in Fig. [5]) is plotted as a function of a for (a) 
PRBUM-COE and (b) PRBUM-CUE. These results are 
fitted with a linear function (solid line) which gives slope 
p = 0.76 and intercept q = 0.72 for the PRBUM-COE 
case. For the PRBUM-CUE case, the slope p = 0.73 and 
the intercept q = 0.77. In both panels, the horizontal line 
represents 7/ = and the vertical line indicates a = 1.0. 



the case of PRBM. However in the PRBM case, one can 
analytically determine the exact transition point by map- 
ping the problem onto a nonlinear a model with nonlocal 
interactions [H] . Because we do not have a good analytical 
tool for treating PRBUM, we need to rely on numerical in- 
vestigations. So we now further investigate the transition 
point from a different perspective. 

For PRBM, the 'not so sharp transition' was attributed 
to the fact that there exists a range of values of a around 
unity in which the eigenstates are fractal up to a char- 
acteristic length denoted by £f rac . This length is anoma- 
lously large as compared to the system size and it diverges 
around a — 1.0 |37j . Following the procedure given in Ref. 
|37j , we now calculate a different kind of fractal dimension 
d.2 of the eigenstates by employing standard box counting 
procedure through the expression 



where S = n/N and X2i n ) is the 2nd moment of the prob- 
ability density of the eigenstate in the boxes of size n. 
In the delocalized (metallic) regime, di should approach 
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Fig. 6. Logarithm of the characteristic length £f rac (de- 
termined from the results in Fig. 5) is plotted as a func- 
tion of a for (a) PRBUM-COE and (b) PRBUM-CUE. At 
am 1.0, log(£f rac ) shows signs of divergence. 



towards 1 from its fractional value as we increase the sys- 
tem size N. So in this regime, the slope of di vs. N curve 
is expected to be positive. On the other hand, in the lo- 
calized regime, the limiting value of c?2 should be 0, and 
consequently the slope of the di vs. N curve will be neg- 
ative. This implies that there exists a critical value of a 
at which the slope will be zero. This particular value of a 
can be regarded as a more precise transition point and on 
that point di will be independent of the system size (i.e. 
scale- invariance), at least for the system size we consid- 
ered. The characteristic length £f rac can be defined as an 
extrapolated system size for which d^ = 1.0 and 0.0 in the 
delocalized and the localized regime, respectively. 

To proceed with this new angle we perform ensem- 
ble averages of d^ over those eigenstates of a matrix that 
have IPR values close to the maximal IPR (thus reducing 
the fluctuations in d-i). Here we construct the ensemble 
by taking 1200 to 300 matrices depends on the Hilbert 
space dimension. In Fig. [5j we plot ^2 as a function of the 
system size N for a = 0.80, 0.90, 0.95, 1.00, 1.05, 1.10 and 
1.20 from top to bottom on a semi- logarithmic scale. The 
results therein can be regarded as a type of finite-size scal- 
ing analysis (quite different from other approaches, e.g., 
[38]). These curves are fitted very well by straight lines. 
We then determine log(£f rac ) as a function of a from these 
fitted straight lines. The results are presented in Fig. [5] 
For both panels in Fig. |6l the left (right) part of log(£f rac ) 
vs a curve shows the fitted results using c?2 = 1 (d% = 0). 
Similar to the PRBM ensemble case [37] , we can see that 
l°g(£frac) diverges around a = 1.0. 

In Fig. [7J we plot the slope r\ determined from the di 
vs. N curves in Fig. 5, for different values of a. The pre- 
sented results suggest that, at a = 1.0, r\ is approximately 
equal to 0.05. This slope is small, but not exactly zero. 
The slope reaches nearest to zero around a = 1.05. In this 
sense, a = 1.05 may be regarded as a more precise tran- 
sition point for both PRBUM-COE and PRBUM-CUE 
cases (applies to b = 1.0 as well). One may suspect that, 
instead of a = 1.0, olq = 1.0 might be the precise tran- 
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Fig. 8. (Color online) Typical eigenvector profile for 
PRBUM-COE [(a)-(c)] and PRBUM-CUE [(d)-(f)]. In (a) 
and (d) a — 0.5; in (b) and (e) a — 1.0, and in (c) and (f) 
a = 1.5. b = 0.1 is fixed. 



sition point. However, according to Fig. [TJ at a — 1.05, 
the parameter ao is approximately equal to 0.95. So we 
have to rule out this possibility. Another interesting ob- 
servation from Fig. 7 is the following: around a = 1.0, 
the slope rj is essentially a linear function of a. The cor- 
responding linear fitting of the form 77 = p — qa gives 
p ~ 0.76 and q ~ 0.72 for the PRBUM-COE case, and 
p ~ 0.73 and q ~ 0.77 for the PRBUM-CUE case. In both 
cases, p ~ q. If we assume p — q in the thermodynamic 
limit, then r) ~ p(l — a), which predicts that rj is zero only 
when a = 1.0. This hints that even here a — 1.05 seems 
to be the more precise critical point for b = 0.1, the slight 
departure from a = 1.0 can be due to the finite size of our 
system. One can also use our PRBUM algorithm to argue 
in favor of a = 1.0 as the true transition point in the ther- 
modynamics limit. That is, if a = 1.0 is not the transition 
point, then our algorithm would have to generate a criti- 
cal PRBUM ensemble starting from a non-critical PRBM 
ensemble. 



5 Statistics of PRBUM eigenvector 
components 

In this section we investigate the statistical properties of 
the eigenvector components of PRBUM. This is interest- 
ing because the parallel results for standard Dyson's cir- 
cular ensembles provide a good reference point for us. 
Figure [5] depicts the profile of some typical eigenvectors 
for PRBUM-COE and PRBUM-CUE, for three different 
values of a, i.e., a = 0.5,1.0, and 1.5, and a fixed b 
value. For a — 0.5 shown in Fig. [5Ja) and Fig. HJd), 
both PRBUM-COE and PRBUM-CUE give delocalized 
random eigenvectors, which is consistent with the obser- 
vation that D2 ~ 1 for a = 0.5. Turning to Fig. [5](b) 
and Fig. HJe) for the case of a = 1.0, the eigenvectors 
exhibit a quite sparse structure, but are not localized in 
any particular regime. The same property was observed in 



the PRBM ensemble [H]. Finally, for the case of a = 1.5 
shown in Fig.[8](c) and Fig. HJf), it is seen that the eigen- 
vectors are highly localized, which is consistent with the 
result of Z?2 ~ seen in Fig. [31 In the following we dis- 
cuss the statistics of the magnitude of PRBUM eigenvec- 
tor components. 

Figure shows the distribution function of the magni- 
tude squared of the eigenvector components for PRBUM- 
COE, considering an ensemble of 100 matrices with N = 
500, for four different a values but with the b value fixed at 
b = 0.1. The results are plotted in terms of the probability 
distribution function P(z), where z = my = ln[iV|(m|</(>A)| 2 ] 
In the case of a = 0.1, because the associated fractal di- 
mensional D2 ~ 1.0 (see Fig.|3J, one expects P(z) to follow 
Dyson's standard COE distribution for eigenvector com- 
ponents. Because the COE distribution of y is given by 2 



p(y) = 



1 



: exp 



(-1). 



the COE distribution for z is given by 

z — exp(z) 



P 



COE 



1 



exp 



(7) 



(8) 



As seen in Fig.^a), the numerical result (open circle) fol- 
lows Eq. © (solid line) remarkably well. The value of a 




Fig. 9. (Color online) Distribution function P(z) vs z 
for the eigenvector component statistics of PRBUM-COE, 
where z — ln[iV|(m|(?i>A)| 2 ] and N is the matrix size. Here 
N = 500 and b = 0.1. (a) a = 0.1, the solid line represents 
the standard COE distribution of Eq. (b) Same as in 
(a) but with a = 0.5. (c) a = 1.0 (around the critical 
point). The central part of P(z) follows closely a normal 
distribution (solid line) with mean value (z) ~ —5.34, and 
variance a z ~ 3.25. Strong deviations from the normal dis- 
tribution are seen for \z\ ^> 0. (d) a = 1.5. An additional 
peak (location indicated by the vertical solid line) appears 
at z = \nN = ln(500) ~ 6.22. This additional peak is due 
to the emergence of some most localized states. 
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Fig. 10. (Color online) Same as in Fig. |H1 but for 
PRBUM-CUE. (a) a = 0.1. The solid line represents the 
standard CUE distribution given in Eq. (fTOj) . (b) Same 
as (a) but for a = 0.5. (c) a = 1.0 (around the critical 
point). P(z) follows closely a normal distribution (solid 
line) with mean (z) ~ —5.93 and variance a z ~ 2.95. 
Strong deviations from the normal distribution are seen 
for z > 0. The left tail (P(z) at z < 0) agrees with the 
normal distribution better than in the PRBUM-COE case, 
(d) a = 1.5. An additional peak (location indicated by the 
vertical solid line) appears at z = In A = ln(500) ~ 6.22. 
This additional peak is due to the emergence of some most 
localized states. 



is then increased to a — 0.5 in Fig. EJb). It is found that 
P{z) is still very close to -Pcoe(<z)j with slight deviation at 
the right tail. This is consistent with our early observation 
that the Di value is still close to unity for a = 0.5. The 
case a = 1.0 is shown in Fig.|9jc), where P(z) is seen to be 
significantly different from the COE distribution Pcoe(z). 
Indeed, the P(z) for PRBUM-COE for this (at least) near- 
critical case is much broader than Pcoe(z), reflecting that 
the possibility of having very large or very small eigenvec- 
tor components is much higher than in the conventional 
random matrix theory. The solid line in Fig. IHIc) repre- 
sents a normal distribution, and it is seen that the central 
part of P(z) is close to the normal distribution. Finally, 
Fig.HJd) displays the result for a — 1.5, a value that yields 
£>2 — > as seen in Fig. [3] In this case, many eigenvectors 
are highly localized and as a result, one observes a much 
enhanced probability at large z. Indeed, because the most 
localized eigenvector has only one component with magni- 
tude unity and all other components with magnitude zero, 
it is now possible to have a high probability at z = hxN. 
The additional peak at z = In 500 ~ 6.22 in Fig. M,d) in- 
dicates that some of the eigenstates are actually close to 
the most localized states. We have also checked that if a 
further increases, then this additional peak becomes more 
and more pronounced. At about a — 10.0 essentially all 
the eigenvectors become the most localized states. 



Let us now perform the same analysis for PRBUM- 
CUE and the results are presented in Fig. [TO] Here the 
reference point is Dyson's CUE statistics of the magnitude 
squared of eigenvector components [2], i.e., 



P(y) = cxp(-y), 



which leads to 



Pcve{z) = exp [z - exp(z)] 



(9) 



(10) 



As seen from Fig. [TOT aL the numerical result for a — 0.1 
agrees with Pcve(z) for the entire range of z. This again 
confirms that PRBUM essentially reduces to CUE for a 
sufficient small a. In the case a = 0.5 shown in Fig. fTOT b). 
P(z) closely follows Pcve{z) with marginal deviation in 
the right tail. The case of a = 1.0 in the delocalization- 
localization transition regime is shown in Fig. HOf cL with 
features clearly different from that of Pcve(z)- We have 
also compared the numerical result with a normal dis- 
tribution [solid line in Fig. ITOT c')]. Interestingly, the left 
part of P(z) as well as its central part is seen to agree 
with the fitted normal distribution reasonably well, but 
clear deviation from the normal distribution is seen for 
z > 0. This feature is different from what is observed in 
the PRBUM-COE case [see Fig. Etc)], where deviations 
from a normal distribution is also clearly seen in the left 
tail. The result for a — 1.5 is shown in Fig. [TOT dL Simi- 
lar to the PRBUM-COE case in Fig. OJd), here again an 
additional peak emerges at z = ln(500) ~ 6.22, which is 
clear evidence that some states become the most localized 
states. 

Summarizing this section, it is observed that when the 
delocalization-localization transition occurs, the distribu- 
tion function P(z) for PRBUM shows strong deviations 
from Dyson's COE/CUE predictions. These deviations 
may be regarded as one main characteristic of the Ander- 
son transition. Interestingly, in the case of PRBUM-CUE, 
both the central part and the left tail of P(z) is close to 
a normal distribution, whereas in the case of PRBUM- 
COE, only the central part of P(z) is close to a normal 
distribution. 



6 Critical spectral behavior of PRBUM 

In this section we investigate the manifestations of the 
delocalization-localization transition of PRBUM in its NNSD 
statistics. In the Anderson tight-binding model on a 3D 
cubic lattice, Braun et al. found that on the metal-insulator 
transition point, the NNSD depends on the boundary con- 
ditions of the cubic lattice [35] ■ Upon averaging over differ- 
ent boundary conditions, it was discovered that the NNSD 
follows approximately the semi-Poisson distribution, i.e., 



P sp (s) = 4sexp(-2s), 



(11) 



where s > is the spacing between two neighboring levels. 
In the s — > limit, P sp {s) ~ 4s, which is a linear level re- 
pulsion property shared by Dyson's COE statistics and 
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Fig. 11. (Color online) The NNSD [denoted P(s)] of the 
eigenphases of PRBUM-COE (open circles) vs level spac- 
ing s, for four different values of a. For the sake of com- 
parison, the Wigner-Dyson distribution for COE (dashed 
line), the semi-Poisson distribution as given in Ea. ([TT|) 
(solid line), and the Poisson distribution (dashed-dotted 
line) are all plotted together, (a) a — 0.1. (b) a = 0.5. (c) 
a = 1.0. (d) a = 1.5. An ensemble of 100 matrices with 
N = 1000 is used for statistics. Results here should be 
connected with those in Fig. [S] 
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Fig. 12. (Color online) Same as in Fig. [TTJ but 
for PRBUM-CUE. Here the dashed line represents the 
Wigner-Dyson distribution for CUE, and the dotted line 
represents the generalized semi-Poisson distribution given 
in Eq. (jnjl with f3 = 2.0. The dashed-dotted line is still 
for the Poisson distribution and and the solid line is still 
for the semi-Poisson distribution. An ensemble of 100 ma- 
trices with N = 1000 is used for statistics. Results here 
should be compared with those in Fig. [TU] 



many time-reversal symmetric dynamical systems with 
chaotic classical limits |1I10) . However, the proportional- 
ity constant in the linear level repulsion for P sp (s) ~ 4s 
is different from Dyson's COE statistics Pcoe(s) ~ f s. 
In the other limit s 3> 0, P S p(s) decays exponentially, i.e., 
P sp (s) ~ exp(— 2s), which is reminiscent of an uncorre- 
lated energy spectrum of a regular quantum system whose 
NNSD follows the Poisson distribution P p (s) — exp(— s). 
Hence in this second limit P S p(s) differs significantly from 
Pcoe(s) which decays to zero as exp(— 7rs 2 /4). Therefore, 
the semi-Poisson distribution P sp (s) can be considered as 
an intermediate critical NNSD behavior. In Ref. [?D], the 
NNSD of the PRBM ensemble was also studied. The find- 
ing is that the central part of the NNSD follows P sp (s) 
well, but non- negligible deviations from P sp (s) were ob- 
served for s — > or for s » regime. We are thus mo- 
tivated to carry out similar NNSD studies for PRBUM 
in the Anderson transition regime and also use the semi- 
Poisson distribution as one of our diagnostic tools. Note 
also that, since Fig. [2] already confirms the uniform distri- 
bution of PRBUM eigenphases, we only need to introduce 
a simple re-scaling factor N/2ir to set the average level 
spacing at unity. 

For the same four PRBUM-COE cases as in Fig. El 
Fig. ITT1 depicts the associated NNSD. In case of a = 0.1 
in Fig. ITTT a). the numerically obtained NNSD (denoted 
by P(s) and represented by open circles) is in good agree- 
ment with the Wigner-Dyson statistics Pcoe(s) from the 
conventional random matrix theory (dashed line). For the 



sake of comparison, in all the four figure panels we have 
plotted the semi-Poisson distribution P sp (s) (solid line) 
as given in Eq. (fTTj) as well as the Poisson distribution 
P p (s) (dashed-dotted line). The case of a = 0.5 in Fig. 
[TIJb) is similar. Consider next the critical case in Fig. 
ITlTc). Quite unexpectedly, P(s) follows the semi-Poisson 
distribution P sp (s) very well in the small-s regime. This 
agreement in the small-s regime is absent in PRBM on the 
critical point. In the large-s regime, deviations from the 
semi-Poisson distribution can be seen clearly in Fig. lllf c). 
Indeed, P(s) in the large-s regime is seen to lie between 
the semi-Poisson distribution and the Poisson distribu- 
tion. This feature is also different from one case study of 
the P(s) of PRBM [40]. In particular, in Ref. [40] the P(s) 
for PRBM actually lies between Pcoe(s) and P sp (s), in- 
stead of between P sp (s) and P p (s). This difference can be 
explained as follows. For b — 1.0 considered in Ref. [4~0l, 
the fractal dimension of PRBM eigenstates is P>2 — 0.66; 
whereas here b = 0.1 and D2 ~ 0.28. As such, the D2 
value in Ref. 10 is closer to unity and the NNSD result 
is expected to be closer to the conventional random ma- 
trix theory. Furthermore, the much smaller D2 value for 
our PRBUM considered here should imply more localized 
states, and hence the P(s) is closer to the Poisson dis- 
tribution. Consistent with this explanation, the case of 
a = 1.5 in Fig. 8(d) has D 2 ~ 0, and the associated P(s) 
indeed agrees well with the Poisson distribution (dashed- 
dot line). 

We next examine the NNSD for PRBUM-CUE in Fig.QU 
The standard Wigner-Dyson distribution for CUE with 
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quadratic level repulsion is now plotted as dashed lines 
in all the four panels. As expected, P(s) in the cases of 
a = 0.1 and a = 0.5 in Fig. HHa) and Fig. [I27b) agree 
with the standard result from random matrix theory. Simi- 
larly, Fig.HHd) shows that P(s) for PRBUM-CUE follows 
the Poisson distribution for a = 1.5, indicating that this 
a value is far beyond the delocalization-localization tran- 
sition point. Much more interesting is the case of a — 1.0. 
Similar to the PRBUM-COE case, in the large-s regime 
P(s) is found to lie between the semi-Poisson distribu- 
tion and the Poisson distribution. For the small-s regime, 
since P(s) in the PRBUM-COE case nicely follows a semi- 
Poisson distribution, one might naively expect that here, 
due to the breaking of time-reversal symmetry, the P(s) 
should follow a generalized semi-Poisson distribution [4"T] 
with quadratic level repulsion. To check if this expectation 
is correct, in all the panels of Fig. we have plotted an 
additional curve (dotted line) to represent a generalized 
semi-Poisson distribution P sp -/3(s) with ft = 2, i.e., 

P sp .^(s)=A^exp[-(P + l)s}, A/3 = ^ 1 - . (12) 

Note that P sp;/ 3 = i(s) = P ap (s). As seen in Fig. lT27 c). in the 
small-s regime, the P(s) for PRBUM-CUE deviates from 
both P sp; ^(s) with [3 = 2 and the semi-Poisson distribu- 
tion P sp (s). This implies the possibility of anomalous level 
repulsion, namely, P(s) ~ s* 3 with (3 being a non-integer. 

In order to investigate carefully the anomalous level 
repulsion for critical PRBUM-CUE in the small-s regime, 
we show in Fig.[T2]the ln[P(s)] vs. ln(s) plot in the small-s 
regime (s < 0.36) and then fit it using the function P(s) ~ 
s* 9 . To show the marked difference between PRBUM-CUE 
and PRBUM-COE around the critical point, this fitting 
is also performed for PRBUM-COE. We obtain /3 ~ 1.04 
for PRBUM-COE, which confirms the linear level repul- 
sion observed earlier in Fig. [TT] By contrast, for PRBUM- 
CUE, we find (3 ~ 1.73, clear evidence of anomalous level 
repulsion. To improve the statistics in the small-s regime, 
we increased the ensemble size by a factor often, i.e., now 
an ensemble of 1000 matrices with TV = 1000 is used for 
our statistics. We have further checked that the change 
in the j3 values thus found numerically is negligible if we 
further increase the ensemble size by a factor of three. 
Thus, at least for the small-s regime we considered, the 
level repulsion is found to be anomalous. In the bottom 
two panels of Fig. [T2J the numerical fitting of the inte- 
grated nearest-neighbor spacing distribution, denoted by 
Z(s), is also presented. The fitting curve I(s) ~ s^ yields 
[3' ~ 2.01 for PRBUM-COE and [3' ~ 2.71 for PRBUM- 
CUE [Note that, by definition of I(s), /3' = [3 + 1]. This 
further supports our earlier observation that P(s) displays 
an anomalous level repulsion in the case of PRBUM-CUE. 



7 Concluding remarks 

The Hermitian PRBM ensemble has been studied exten- 
sively in the literature as a standard model for the An- 
derson transition. In terms of the power-law decay of the 




ln(s) ln(s) 



Fig. 13. (Color online) The NNSD [denoted P(s)] and 
the integrated nearest neighbor spacing statistics [de- 
noted I(s)] in the small-s regime for both PRBUM-COE 
and PRBUM-CUE. In the top two panels, the numeri- 
cal data (open circles) are fitted by P(s) ~ s&. In the 
case of PRBUM-COE, [3 ~ 1.04 (almost linear repulsion), 
whereas in the case of PRBUM-CUE, [3 ~ 1.73 (signature 
of anomalous level repulsion). For the bottom two pan- 
els, I(s) vs s is fitted by I(s) ~ s' 3 ', with ~ 2.01 for 
PRBUM-COE and /?' ~ 2.71 for PRBUM-CUE, thereby 
further confirming the found j3 values. To improve the 
statistics in the small-s regime, an ensemble of 1000 ma- 
trices with N = 1000 is used for our statistics. 



magnitude of the matrix elements from the diagonal, the 
unitary PRBUM ensemble is similar to PRBM (the gen- 
eration of PRBUM is much more involved). This work 
showed that, statistics of PRBUM in the delocalization- 
localization transition regime can be different from their 
PRBM counterparts. Hence PRBUM is not a trivial uni- 
tary analog of the Hermitian PRBM ensemble. Somewhat 
consistent with this conclusion, we have observed that for 
a given set of b and a, the matrix elements of PRBUM as 
generated from our algorithm decay from the diagonal as a 
different power-law function [see Eq. @ and Fig. 1, where 
bo and ao can be appreciably different from b and a] . Given 
our detailed findings, we hope that PRBUM can capture 
some critical behavior of actual Floquet eigenstates in the 
delocalization-localization transition regime, by modeling 
system's periodic unitary evolution matrix with power- 
law decaying matrix elements. Different choices of ensem- 
ble parameters can model a system in either insulating or 
metallic phase, and if the ensemble parameters are close to 
critical values, then fractal behavior of the Floquet eigen- 
states (within a certain scale) may be also modeled by 
PRBUM. One obvious advantage of PRBUM as compared 
with PRBM is that in performing the statistics all the 
eigenstates are treated under equal footing, rather than 
selecting certain energy windows as often done in PRBM 
studies. We also note that in our early study [27 PRBUM 



10 



Jayendra N. Bandyopadhyay, Jiangbin Gong: Title Suppressed Due to Excessive Length 



can correctly model how the variance of D2 of Floquct 
eigenstates (in a double-kicked top model) scales with the 
system size but PRBM cannot. 

In summary, we have numerically studied the statisti- 
cal properties of PRBUM in the derealization- localization 
regime and have compared them with PRBM and the 
standard random matrix theory when feasible. The de- 
pendence of the fractal dimension of PRBUM eigenvectors 
upon ensemble parameters is used to examine where the 
delocalization-localization transition point is, within some 
uncertainty due to finite-size effects. For cases very close 
to the transition point, the central part of the PRBUM 
eigenvector-component distribution function is close to a 
log-normal distribution and the tails (especially the PRBUM- 
COE case) show strong deviations from the log-normal 
distribution. We have also examined the NNSD of PRBUM 
and have compared the result with the semi-Poission dis- 
tribution and a generalized semi-Poisson distribution. An 
anomalous level repulsion of PRBUM in the absence of 
time-reversal symmetry is also numerically found. 

This work was supported by the Academic Research 
Fund Tier I, Ministry of Education, Singapore (grant No. 
R-144-000-276-112). 
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